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ABSTRACT 


A method is described which was used to measure mu- 
on energy spectrum characteristics in muon groups 
underground using q-e decays recording. The Baksan 
Telescope’* s experimental data on p-e decays intensi- 
ty in muon groups of various multiplicities are 
analysed. The experimental data indicating very flat 
spectrum does not however represent the total spect- 
rum in muon groups. Obviously the muon energy spect- 
rum depends strongly on a distance from the group 
axis. The "core attraction" effect makes a significant 
distortion , making the spectrum flatter .After taking 
this into account and making corrections for this ef- 
fect the integral total spectrum index in groups has 
a very small dependence on muon multiplicity and ag- 
rees well with expected one: 6 = e e xpected = :IL * 75 • 


1 . Introduction. It has been shown (1) that p-e decays method 
can be used to measure an effective exponent of muon energy 
spectrum in muon groups underground, which is expected (2) 
to be much flatter than the total muon flux spectrum*). 
Unfortunately, such measurements are affected by so called 
"core attraction" effect, which preferably selects the cen- 
tral part of the muon group making the spectrum of recorded 
muos flatter, as compared with the spectrum in the group as 
a whole. This effect depends on the ratio of detecting area 
size to the muon group dimension, on the form of lateral dis- 
tribution function and on muon multiplicity spectrum. To make 
the corrections for "core attraction" effect we have to use 
the above mentioned, parameters which we take from (1) and 
(7). Because the exponent of multiplicity spectrum is not a 
constant at small multiplicities , we have to do the analysis 
for different multiplicities separately. 

For the total muon flux the relations between experimen- 
tally observed quantities R a ^ and Ri oc and muon energy spec- 
trum underground are well established and understood (3,4,5). 
By definition, R a j. is the ratio of stopping to throughgoing 
muon fluxes normalized to 100 g/sm 2 target. R lpc is defined 
in similar way, but for cases when muon from atmosphere does 


*) First experimental evidence of this has been presented 

in (5) . 
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not stop and p-e decays come from pions originated in muon 
induced hadronic cascade. 

/!/; = /00. a -./ n'/ 

^ ( >0 ) ° -OP 

<’<£/• w > nn **• 13 o-e?”(S-o-rJ ’ /2V 

d3(s„)/JSf, * (s^sj,)^' 1 m 


The eqs. 11-21 give these relations in general way. Here, 
£*, is muon energy underground, a=2.34>10 _ 3 Gev*sm 2 /g is mean 
ionization energy losses at our depth, c=2*10“5 Gev~-75 (ad- 
justed by experiment ). The eqs. /I*/ and / 2 * / are valid if 
the spectrum underground is expressed by eq. / 3 / , £ 0 is a 
specific energy equal to 210 Gev for our depth. Eq. / 3 / 
corresponds to purely power law muon spectrum on the 

surface. Nevertheless, if the spectrum is not strictly a power 
law eqs. /I’/ and /2'/ can be used to determine the effec- 
tive exponent y near the depth energy threshold. But, in the 
case of muon groups such a procedure would be erroneous be- 
cause of distortion due to "core attraction". 


2 . Experiment . Iflnthis experiment the u-e decays were recorded 
in the second horizontal 200 m 2 and 40 tons scintillation la- 
yer of Baksan Telescope (4). Other 7 layers were used to de- 
termine the number of parallel muon trajectories (m) and to 
distinguich between "atmospheric" and "locally produced" 
events. The latter procedure is made with high confidence 
(~99%) for low multiplicities , but becomes difficult for m>30. 
The selection of m> 1 events was made by the off-line compu- 
ter. All necessary data associated with these events were 
printed and analysed visually together with 10-beam oscillo- 
scope slides of decay's electrons pulses. The only trigger 
used in the experiment was a delayed pulse in the mentioned 
scintillator layer with rate «30 h'^.Data accumulated during 
200 h of observation were included in the analysis. A flux of 
muons throughgoing the target layer in groups has been mea- 
sured in a special short run with the same conditions and da- 
ta processing but without requirement of p-e decay. 

3. Calculation of R a f-. and R -| no in muon groups. To do this the 
assumptions as in (1,2,73 and a number of consecutive integ- 
rations have to be made: 


f(A/) tfV / 4 / ; A/i //• 4 (i) / 5 / ; a(H ) = jj'ffrj J$j jfrz F(z)F Z / 6 / 

/r p \0.1 °o 

F *(E h t) - Ey*e & "£) /7/; F(*)^ [F‘(£/,,i)etE/* W 


. elk! 




/9/ 
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Here in / 4 / : N is the mean number of muons at a given 
depth from definite primary, y(N)dN is a distribution of N 
assuming standard energy spectrum and chemical composition 
of primaries. 


In / 5 / : A(R) is a fraction of N which happened to be 

inside detecting area,R is a distance between centre of area 
and muon group core. 

In /6 / : A(R) is calculated using lateral distribution 

function F(r), where r is a distance between area element dS 
and core position. 

Eq. /7/represents the assumed lateral-energy distribution 
of muons, E^ is muon energy at the surface, 6 is integral ex- 
ponent of power law muon spectrum there, E t is the minimal 
energy to penetrate to a given depth. For our depth r n =6.2 m 
if 6 = 1 . 7 5 . Eq . /8/ represents muon lateral distribution. 

In /9/:Gjp(R) is distribution of impact parameter R for a 
given multiplicity m. Poisson fluctuations of m for a given N* 
are assumed. Fig. 1 shows these normalized distributions for 
several m. Curve (m^l) represents the case without core att- 
raction: G^i(R)~ 27T*R*A(R) .One can see that for m>3 the dis- 
tortion of distribution of R becomes standard, for m=l the 
distortion is small, but of the opposite sign. 

In /10/:by integrating F»(E M ,r) over detecting area S, r 
being the distance between dS and the group core, then integ- 
rating over R we get the muon energy spectrum inside the de- 
tecting area for a given recorded multiplicity m. 

Eq. /II/ represents the normalized muon energy spectrum 
at a given R.Eq. /12/ is relation between muon energy under- 
ground (£y) and on the surface (E w ). 

Next steps are as follows :using /12/ transform /10/ into 


spectrum underground ;put this spectrum into eqs./l-2/ and ob- 
tain the expected quantities R„ t and R loc for assumed 8; se- 
lect 6 which fits the experimental values of R . and R lo .So, 
for each m we obtain two values of 8: 6^ related to R afc ,and 
62 related to R-^ 0 .If all the procedure is reliable and the 
surface muon spectrum follows power law, then 6-iand 6 P should 
coincide. 


4. Results. All the data are presented in Table. Denoted as 
* * 

6^ and indices are exponents obtained without taking into 
account the "core attraction" effect, but simply using eqs. 
/ 1 *t an d /2*/. It is seen^that corrected exponents 8 are 
bigger than noncorrected 6 . On the other hand, one can see 
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DISTANCE FROM DETECTOR CENTRE . R/RB CRB-fi.2 M5 

Fig.l. Distribution of inpact parameter R for various multiplicities. 

that 6 . and 6 2 are closer to each o.ther than 6 . and 6 2 which 
indicates a good fit to the power law spectrum, what can be 
expected. Deviations from this for m=ll+30 has probably sta- 
tistical origin, for m =2 it could be expected because of pri- 
mary nucleon energy for m =2 is not big enough (7). 

Taking all events with m^-3 and using both 8 . and 6 2 we 
obtain the integral spectral index of muons in groups as 
8=1.67+0.11 for E > 200 Gev . This result is in agreement 
with EAS data ( 6 ). M 

Table 


m 

1 2 

3+4 

f 5+10 

n +30 

I>1 

>2 

>3 

R at x103 

■ 

gQjjjjfl 

1.29 
+ .26 

■ 

1 

0.66 

+.29 

| 

I 

1.17 
+ .09 

1.16 
+ .17 

R loc xl0 ‘ 

■ 

1.27 
+ .23 

1.75 
+ . 38 

~T7&B 

+ .47 

0.75 
+ .03 

1.27 
+ • 08 

1.76 
+ .20 

* 

6 1 

B 

1.17 
+ .24 

1.05 

+.30 

0.60 
+ . 26 

'2.10 
+ .08 

1.06 

+ •08 

■ 

1.05 
+ .15 

* 

82 

2.08 

±.25 

1.56 
+ .20 

1.27 
+. 16 

1.31 
+ .23 

2.36 
+ .08 


1.27 
+ .08 

6 1 

1.55 
+ .14 

1.83 
+ .27 

1.73 

+.35 

1.18 
+ .41 

2.10 
+ . 08 


1.71 
+ .25 

&2 

2.32 
+ .21 

B9| 

s 


1.73 
+ . 18 

2.36 
+ .08 

I 

1.66 
+ .13 
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